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Abstract 

We show one can use classical fields to modify a quantum optics ex- 
periment so that Bell's inequalities will be violated. This happens with 
continuous random variables that are local, but we need to use the cor- 
relation matrix to prove there can be no joint probability distribution of 
the observables. 
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1 Introduction 

The issue of the existence of hidden variables for quantum mechanics is almost as 
old as quantum mechanics itself. However, in 1963 J. S. Bell showed j^, ^| that 
if one makes some "reasonable" assumptions about the hidden variables, like 
locality and statistical independence of distant measurements, the correlations 
for the outcome of measurements for an EPR-like experiment have to satisfy a 
set of inequalities. In 1982 Alain Aspect and coworkers showed that quantum 
mechanics violated Bell's inequalities, drawing the conclusion that one cannot 
have a local realistic theory that would replace quantum mechanics |^|, 

Because Bell's assumptions were considered equivalent to the existence of 
an underlying physical reality, it is often said that any classical system satisfies 
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Bell's inequalities. In this paper we will show that classical fields do not satisfy 
Bell's inequalities, hence classical fields, e.g. electromagnetic fields, are not Bell- 
type hidden variables. We do this by showing that a simple experimental setup, 
suggested by Tan et al. |], ^| , can be reinterpreted for classical electromagnetic 
fields. For this reinterpretation we derive from the classical field properties a 
violation of Bell's inequalities H f|, |tJ, with, at the same time, locality being 
preserved in a sense to be made precise. 

2 Experimental Setup 

The experimental scheme uses two classical coherent sources a\{6\), with phase 
6\, and 0^2(^2), with phase #2, and a third source to be studied, u(9), with 
unknown phase. The experimental configuration has two homodyne detections, 
(Di,D2) being one and {D^,Dn) the other, such that the measurements are 
sensitive to phase changes in u(9). The geometry of the setup is shown in FIG. 
1. In FIG. 1 BSl, BS2 and BS3 are beam splitter mirrors that will reflect 
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Figure 1: Proposed experimental configuration. 

50% of the incident electromagnetic field and let 50% of it pass. When the 
electromagnetic field is reflected, the mirrors add a phase of 7r/2 to the field, 
while no phase is added when the field passes through BSl, BS2 or BS3. We 
will look for correlations between the pairs of detectors (D\, D2) and (D3, D4). 

3 Correlation Functions 

In this section we compute the correlation functions that violate the inequalities. 
We first define the random variables in terms of which we derive the Bell-type 
correlations. On this matter we shall be as explicit as possible. Associated to 
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the source u(9) at D\ is the random variable Ui(t), whose value at t is just the 
value of the classical field at D\, namely, 

U 1 (t) = -j3cos(ujt + 6 + n/2), (1) 

where (3 is the amplitude of the field at the source, 9 is the unknown phase and 
7r/2 is a phase gained when u is reflected at BS3. 

Probability enters by using the time average to compute the expectation of 
Ui(t) 2 

U 2 = (U^tf) = ([i/3cos(^ + + 7r/2)] 2 ), (2) 

which is just the standard intensity, but here we treat it probabilistically. In 
a similar fashion, associated to the source a\{9\) at D\ is the random variable 
Ai(t), 

A 1 (t) = ^O!cos(wt + 0i +7T/2) (3) 

and thus 

A\ = (Ai(t) 2 ) = ([±acos(wi + 9 1 + n/2)} 2 ). (4) 

At Di, the total field is the random variable F\(t) = U\{t) + A\(t). So, the 
intensity of the total field at D\ is just the second moment of F\(t), i.e., 

h(9) = F 2 = {F l (tf) = {{U l {t)+A l {t)f) 

= (U 1 (tf)+2((U 1 (t)A 1 (t))) + (A 1 (tr), (5) 

where we used 9 as an argument for I\ to make it explicit that it depends on 
9. We can see that the cross moment in the expression above is the classical 
interference term. 

We can compute I\ directly from the expression for Ui(t) and A\(t) in the 
following way 

h{9) = lim i / [Jacos(wt + 0i+7r/2) + 

T-»oo 1 J Q I 

i/3cos(wt + 6» + 7r/2)] 2 dt, (6) 

which is 

h{9) = ^ 2 + \^cos{9-9 l )+ 1 -a 2 . (7) 
In similar fashion, we can compute for the other three detectors, 

h{9) = ^p 2 -\ a pcos{9-9 1 )+ l - a 2 , (8) 
h{9) = ^l3 2 -\ap S \n{9-9 2 )+ l -a 2 1 (9) 
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and 



h{0) = + ~a/3 sin(0 - 9 2 ) + ±a 2 . 



(10) 



The intensities obtained are conditional on 9. To obtain the unconditional 
intensities we assume a uniform distribution for 9 and integrate the expressions 
for all possible values of 9. Not only is 9 unknown, but the phase would vary 
randomly for repeated runs of the experiment. If 9 were a coherent source with 
fixed 9, Bell's inequalities would not be violated (§]. 

The unconditional intensities Ii, I2, I3, and I4 for the detectors D\, D 2 , -D3, 
and D4 are 



h = 
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We can see from that the intensities are the same for all detectors, 

and are similar to those given by Walls and Milburn 0] in the case of a classical 
source. 

We now start computing the covariance between intensities in the homodyne 
detectors. The covariance we are interested in is between (I\ — I 2 ) and (^3 — I4). 



Cov(/ 1 



h,h-h 



1 

27 



[{h{9) - I 2 {9)) x (I 3 (0) - h{6))]dS 



- ^ J \h(?) - h{6))d6 x i- jf VsW - 



It is straightforward to show from (|7|)-(|l0|) and (|l5| ) that 
Gov{h -I 2 ,h- h) = -^/? 2 « 2 sin(0i - 



2). 



(15) 



(16) 



In order to compute the correlation we have to know the variance of the 
random variables (7i — 1%) and (I3 — I4), which are defined as 



Var(Ji-I 2 ) = ^ JJ {h{6)-h{6)Yd9- i- J Q \h{9) - h{6))d6 



32 M 
(17) 



and 



Var(7 3 -/ 4 ) 



2tt 



(J 3 (0)-/4W) 2 ^- 



1 

2^ 



2tt -|2 

{h{9) - h{d))d6 



II 
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(18) 
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Finally, we are in a position to compute the correlation between the two random 
variables (I1—I2) and (I3— 14). This is done in the standard way, by just dividing 
the covariance by the squareroot of the variances: 

P {h hJs - h) = Covjh- I 2 ,I 3 -h) 

V Var(/! - I 2 ) Var(/ 3 - h) 

and we have the following expression for the correlation 

p{h - h,h ~h) = - sm(0 x - 2 ), (20) 

which we may rewrite as 

p(e 1 ,e 2 ) = -sia(9 1 -9 2 ). (21) 



4 Violation of Bell's Inequalities 

We are now in a position to show that we can violate Bell's inequalities. We 
may now choose angles 6\, 62, 0[, and 0' 2 such that we obtain at once, for the 
four correlations p(0i,02), p(0i,0' 2 ), p(0' 1 ,0 2 ) and p(0' 1 ,0' 2 ) a violation of Bell's 
inequalities in the form due to Clauser, Home and Shimony by choosing the 
four angles such that 

9i-9 2 = 0[ -9' 2 = 60°, (22) 
0i- 02 = 30°, (23) 
0[ - 2 = 90°. (24) 



In particular, 



p(9i,9 2 ) - P i0iA) + p(0'M + Pie'iA) = 

-^ + 1-1-^ < 2. (25) 
2 2 2 y ' 

However, in the case of continuous random variables, which is what we have 
in the present context for intensity, or differences of intensity, failure to satisfy 
Bell's inequalities in the Clauser, Home and Shimony form does not imply that 
there can be no joint distribution of the four random variables compatible with 
the four given correlations. In fact, it is easy to show that for selected values of 
the two missing correlations, there does, for this example, exist a joint proba- 
bility of the four random variables compatible with the four given correlations. 
What is required in the general case, as opposed to that of discrete ±l-values, 
to test for the existence of a joint distribution, when means and correlations 
are given, is that the eigenvalues of the correlation matrix are all nonnegative. 
This extends the earlier result of for discrete ±1 values. Because of the 
freedom to select arbitrarily the two missing correlations in the Clauser, Home 
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and Shimony form of Bell's inequalities, we have not been able to construct an 
example using — sin(0, — 6j) for the correlations that has at least one negative 
eigenvalue for all possible values of the missing correlations. 

Another possibility is obvious. Bell's original paper || used three rather 
than four random variables, and, put in terms of this letter, he showed that the 
correlations for the three angle differences derived from 1; 9 2 , and 9 3 violated 
the following inequality, necessary for the existence of a joint distribution of 
three discrete random variables with values ±1: 

P (9 1 ,9 2 ) + p(0i,6a) + P (9 2 ,9 3 ) > -1. (26) 

To violate ( p6| ) we choose three angles 9%, 9 2 and 6*3, with 

9i = 0, (27) 

2 = 45°, (28) 
9 3 = 90°, (29) 
and, using equation (|Tj), the correlation matrix is 

1 — y/2/2 -1 \ 
-a/2/2 1 -y/2/2 ■ (30) 
-1 -y/2/2 1 / 

It is a direct computation to show this matrix has both positive and neg- 
ative eigenvalues, so that it is not nonnegative definite. The eigenvalues are 
(y/b + l)/2, (-V5+ l)/2 and 2. Therefore, there can be no joint probabil- 
ity distribution for the three random variables compatible with the correlations 



given in (30). From results in |llj] and [12J there can then be no hidden variable 
that factors out the correlations conditionally, i.e., there can be no A such that 
for the three random variables X(9i), Y(9 2 ) and Z(9 3 ), we have 

E{XYZ\X) = E(X\X)E(Y\X)E(Z\X), (31) 

since there is no joint probability distribution of X , Y and Z compatible with 
the given correlations. In particular, 9 = A cannot serve as a Bell-type hidden 
variable for a classical field described by (3). 

Finally, we note that even though (E6) was violated by the angle values in 
(|7|)-(|9|), this inequality is not a satisfactory general test for existence of a joint 
distribution, as the following shows. Let 9\ = 0, 9 2 = 30° and #3 = 45°. Then 
it is easy to check that inequality (^) is violated, but the eigenvalues of the 
correlation matrix are all nonnegative, and so a joint distribution exists. 
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5 Measurement and Photon Counts 



Because classical field theory is a deterministic theory, our introduction of expec- 
tations and probabilities might be questioned. Our response is that the strength 
of a classical field at a space-time point cannot be measured, as was emphasized 
long ago by Bohr and Rosenfeld in a famous paper in 1933 Jll|. As they pointed 
out, classical field strength cannot be represented by true point functions, but 
by average values over space-time regions. This is exactly what we have done in 
introducing random variables and their expectations. The casual reader might 
claim that we should do an analysis of coincidence counts with photocounters. 
This makes no sense in the case of classical fields, where the number of photons 
arriving at the same time at each detector is incredibly large. What makes sense 
is not discrete but continuous measurement of intensity. 

Despite that, we are going to use the previous result to model discrete photon 
counts in such a way that they violate Bell's inequalities. For this, we define two 
new discrete random variables X = ±1 and Y = ±1. These random variables 
correspond to nearly simultaneous correlated counts at the detectors, and are 
defined in the following way. 

„ _ J +1 if detector D\ triggers a count , . 

\ — 1 if detector Di triggers a count ^ ' 

y _ ( +1 if detector D 3 triggers a count , , 

\ — 1 if detector D4 triggers a count. 

To compute the expectation of X and Y we use the stationarity of the process 
and do the following. First, let us note that 

h-h = N x -P{X = l)-N x -P{X = -\), (34) 

where Nx is the expected total number of photon counts at D\ and D2 and 
P(X = ±1) is the probability that the random variable X has values ±1. The 
same relation holds for 

I 3 -I 4 = N Y -P(Y = l)-NyP(Y = ~l). (35) 

To simplify we put as a symmetry condition that Nx — Ny = N, i.e., the 
expected number of photon counts at each homodyne detector is the same. But 
we know that 

h+I 2 =N-P(X = l)+N-P{X = -l) = N, (36) 

and 

I 3 + h = N ■ P(X = l) + N ■ P(X = -1) = N. (37) 

Then we can conclude from equations (0) — (10), assuming maximum visi- 
bility, that 

E d (X\8)= I -±—^=cos(9-9 t ), (38) 
Ji + H 
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E d (Y\e) = 1 -^— ^ = 8^(9-6^, (39) 

where E d represents the expected value of the counting random variable. It is 
clear that if 9 is uniformly distributed we have at once: 

E{X) = E 9 (E d {X\6)) = 0, (40) 

E(Y) = E B {E d {X\0)) = 0. (41) 

We can now compute Cov(X,Y). Note that 

Cov(X,Y~) = E{XY) - E(X)E(Y) 

= E e {E d {XY\9)) - E e {E d {X\9))E e {E d {Y\9)) (42) 



Cov(X,r) = — / E d (XY\9)d9 





1 i-27T 



and so 



, . E d {X\6)d6 x — / E d {Y\9)d6. (43) 

In order to compute the covariance, we also use the conditional independence 
of X and Y given 9 7 which is our locality condition: 

E d (XY\9) = E d (X\9)E d (Y\6), (44) 

because given 0, the expectation of X depends only on 9i, and of Y only on 6j. 
Then, it is easy to see that 

p(X, Y) = Cov(X, Y) = — sin(0* -6j). (45) 

The correlation equals the covariance, since X and Y are discrete ±1 random 
variables with zero mean, as shown in (|4^) and (|4l|), and so Va,r(X) = Var(Y") = 
1. It follows at once from ( |45| ) that for a given set of 0j's and Oj's Bell's inequal- 
ities are violated. 



6 Locality 

Much of the discussion involving Bell's theorem is connected to locality. For that 
reason, we will prove in this section that the scheme presented in this paper is 
local in one precise sense. We follow j7j. Locality requires the following: 

E(X\9 1 ,9 2 ,6) = E(X\9 ll 9). (46) 

It is obvious that (|4^) follows immediately from subtracting (||) from ([?]), and 
observing the result does not depend on #2, and similarly for the other cases 
of random variables Y and Z. Equation ([hi|) says simply that whatever is the 
result of the measurement at one homodyne detector, it must depend only on 
9, the hidden variable, and the phase associated to this particular detector, and 
cannot be influenced by the phase at the other detector. 
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7 Proposed experiment. 



The experiment proposed in pj] supposes a single photon source that is split into 
the two homodyne detectors. Tan et al. also analyze the classical case and get 
no violation of Bell's inequalities. However, they assume a weak coherent source 
with randomized phase as the classical analogue of their single photon source. 
This would be equivalent to having a classical thermal source, where coherence 
would not be a strong feature. In our experiment we suppose that this source 
is not only classical, i.e., with high intensity, but also that it is coherent with 
the phase unobservable and varying randomly on repeated runs. The different 
source used here, as opposed to that used in [jj] implies that the expectations 
given by (7)— (14) are computed in a different way than in [jj]. Here we first 
integrate with respect to t and them integrate with respect to 9. It is easy to 
supply a source that would fit our requirements. This would be, for example, 
a radio source, a microwave, or a laser source, all with unstabilizcd phases. 
To realize this experiment, one must also use two additional coherent sources 
with stable known phases and with the same frequency as the nonstabilizcd 
source. If a data table is then built that keeps track of all the measured values 
on the detectors, we can compute the correlations and see a violation of Bell's 
inequalities, or do the stronger test using the 3-variable version and the matrix 




8 Final Remarks 

There are several remarks that we must add in order to clarify some points. 

First, when using classical fields the number of photons is overwhelmingly 
large. For that reason, we would not need to compute any photon count corre- 
lation. What we measure is intensity. On the other hand, Bell's inequalities are 
not enough to show that we do not have a joint probability distribution for clas- 
sical fields, because they assume a continuous range of values. That is why we 
computed the correlation matrix, showing that for this case a joint probability 
distribution does indeed not exist. 

Another point is that intensity of classical fields does not satisfy the basic 
assumption made by Bell, because it can take an infinite range of values; Bell 
considered spin measurements that can take only two possible values. To show 
that this does not present any constraint, in Section V we did an analysis of 
photon counts, which can only take, as in Bell's assumptions, two discrete values. 

Finally, the last point. One can argue that if classical fields violate Bell's 
inequalities, then, since they are classical, Bell's theorem must be wrong, and 
we must show why it is wrong. We did not show that Bell's theorem is wrong. 
We just showed that a classical field is not a Bell-type hidden- variable. What is 
wrong is the preconception that anything classical must satisfy Bell's hypothesis. 
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